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ABSTRACT 

This paper presents a detailed study of the scattering polarization profiles 
formed under partial frequency redistribution (PRD) in two thermal models of 
the solar atmosphere. Particular attention is given to understanding the influ- 
ence of several atmospheric parameters on the emergent fractional linear polar- 
ization profiles. The shapes of these Q/I profiles are interpreted in terms of 
the anisotropy of the radiation field, which in turn depends on the source func- 
tion gradient that sets the angular variation of the specific intensity. We define 
a suitable frequency integrated anisotropy factor for PRD that can be directly 
related to the emergent linear polarization. We show that complete frequency 
redistribution is a good approximation to model weak resonance lines. We also 
show that the emergent linear polarization profiles can be very sensitive to the 
thermal structure of the solar atmosphere and, in particular, to spatial variations 
of the damping parameter. 

Subject headings: line : profiles - polarization - radiative transfer - scattering - 
Sun: atmosphere 
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Introduction 



The linearly polarized sola r limb spectrum produ ced by scattering processes in quiet 
regions of the solar atmosphere (jStenflo &: Kellerlll997l ). contains a wealth of information on 
the physics of scattering and on the thermodynamical and magnetic conditions of the solar at- 
mosphere. The rich structuring of this so-called second solar spectrum and its i nterpretation 
has opened a new window in solar phy sics with great diagnostic potential (e.g.. IStenflol 12004 
20061 ; iTrujillo Buenoll200ll . 120031 120091 ) . This linearly polarized spectrum h as been measured 
with high s pectral resolution from the UV at 3160A to th e red at 6995 A (Gandorfer 2000. 



2002 



20051 ). Using this atlas of the second solar spectrum iBelluzzi fc Landi Degl'Innocenti 



( 120091 ) have presented a very useful classification and spectroscopic analysis of the most 
polarizing atomic lines. 

The modeling and physical interpretation of the second solar spectrum requires the 
application of sophisticated theories of line formation capable of accounting for, in a self- 
consistent way, the physics of scattering polarization in the presence of magnetic fields. 
Starting from the principles of quantu m electrodynamics, a rigorous theo ry of the polariza- 
tion in spect ral lines was formulated by iBommier fc Sahal-Brechotl (119781 ) for optically thin 



lines and by lLandi Degl'Innocentil (119831 see also the monograph by Landi Degl'Innocenti 
& Landolfi 2004) for optically thick lines. This latter theory, which can take into account 
radiative transfer effects in multi-level atomic systems, lower level polarization, hyperfine 
structure and level crossings interferences, has be en successfully applie d to interpret several 
interesting spectropolarimetric observations (see ITrujillo Buenol 120091 . for a recent review). 
However, a limitation of this theory is that it is based on the approximation of complete 
frequency redistribution (CRD) in scattering (namely, no correlation between the frequen- 
cies of the incident and scattere d photon). Though this is a suitable approximation for 



several solar spectral lines (e.g.. iManso Sainz fc Trujillo Buenol l2003t ITrujillo Bueno et al. 



2004 ; IStepan fc Trujillo Buenol 120101 ). it is not suitable to model the observed linear po- 
larization patterns in some of the strongest and most prominent lines of the second solar 
spectrum, like Cal 4227 A and Call K at 3933 A. For these and a few other resonance 
lines one has to take into account the effects of partial frequency redistribution (PRD) in 
scattering. 

An attempt to handle PRD effects within a theoretical framework where the basic 
results of the CRD density-matrix theory are generalized through the hypothesis of consid- 
ering the atomic levels as a continuous dis tribution of infinitely sharp s ubleye ls (the so-called 



metalevels theory) has been proposed by lLandi Degl'Innocenti et al.l ( 19971). Although this 



theor y has proved to be quite successful for some applications (e.g., 



Landi Degl'Innocenti 



19981 ) the difficulty of coherently incorporating depolarizing collisional effects in such a theo- 
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retical approach prevents us from employing it in the present investigation. We thus prefer, 
in this paper, to rely on the better-established theory of resonance scattering (see below) 
although we have to remind that such a theory suffers from the main limitation of being 
applicable only for two-level atoms with unpolarized, infinitely sharp lower-levels. 

In a scattering event, the direction, frequency, and polarization of the scattered photon 
are in general different from those of the incident photon. The correlations between angle, 
frequency and polarization of the incident and scattered photons can be described by a PRD 
matrix. The effects of PRD are clearl y observable in the wings of strong resonance lines and 
to a lesser extent in the line core (e.g.. lFrischlll996l ). The general proble m of redistribution o f 



resonance radiation including the effects of collisions was investigated by lOmont et al.l ( 119721 ). 
Using a quantum mechanical description of matter and radiation, they derived PRD functions 
in the rest frame of the atom. A year later, t hese authors addressed the s ame problem 



but fo r the magnetized case (|Omont et al.lll973l ). Starting from the work of lOmont et al. 



( 119721 ). iDomke fc Hubenyl (119881 ) derived expressions of the PRD matrices for resonance line 
polariz ation in a two-level atom with unpolarized lower-level. By applying a master equation 
theory, iBommierl ( 1l997a[ ) derived a more elegant but equivalent expression for the non- 
magnetic PRD matrix. Furthermore, she derived the explicit form of the labo ratory frame 



PRD matrix in the presence of an arbitrary magnetic field (see iBommier 



1997b). for the case 



of a tw o-level model atom without atomic polarization in the lower level. IBommier fe Stenflo 
( 119991 ) developed a classical time-dependent oscillator PRD theory for the particular case 
of a normal Zeeman triplet in the atomic rest fr ame. The ensuing labora tory frame PRD 
matrices for the magnetized case were derived by ISampoorna et al.l (l2007al Jbl). In this paper 
we restrict our attention to the problem of resonance line polarization in situations where the 
radiation field is axially symmetric (e.g., the case of non-magnetic resonance line polarization 
in plane-parallel stellar at mospheres). We use the angle-averaged (A A) version of the non- 
magnetic PRD matrix (see lDomke fc Hubenyl Il988l ; lBommierlll997bl ). 



Dumont et al.l ( 119771 ) carried out model calculations of polarized line transfer with PRD 



in axi-symmetric plane-parallel atm ospheres. They used the angle-dependent type I (pure 
Doppler) redistribution function of iHummerl (119621 ). which may be used in some circum- 
stances for the Doppler cores of resonance lines. They considered (1) an isothermal model 
atmosphere and (2) a chromospheric-like model atmosphere described by an analytic form 
for the Planck function. For the case of isothermal model atmospheres they studied the 
influence of e (the photon destruction probability per scattering event) and the error intro- 
duced by the assumption of CRD on the emergent linear polarization profiles. They showed 
that CRD is adequate to describe the line core polarization. For the case of schematic 
chromospheric-like model atmosph eres, they cons i dered the influence of the line strength, 
but under the assumption of CRD. iRees fc Salibal (119821 see also Saliba 1985, 1986) consid- 
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ered t he same problem, but using an AA type II redistribution func tion (-R t t.aa) of lHummer 
( 119621 ). However, for -Ri^aa they used an approximate form given by iKneerl (ll975[ ). which as- 
sumes that CRD prevails in the line core and coherent scattering in the wings. They studied 
the influence of the line strength, e and damping parameter a of the absorption profile on the 
emergent linear polarization and showed that in the presence of a background continuum, 
linear polarization profiles with both core and wing maxima are formed in isothermal as well 
as in chromospheric-like model atmospheres. 



Later, iFaurobertl ( 119871 . Il988l ) studied the linear polarization of resonance lines formed in 
isothermal slabs of finite thickness, in semi-infinite isothermal atmospheres and in chromospheric- 
like model atmosphere. She considered both angle-dependent and AA type II redis tribution 
functions of iHummerl (ll962f ). She confirmed the conclusion of iDumont et al.l ( 119771 ). namely, 
that the approximation of CRD is adequate to describe the line core polari zation even when 
type II redistribution is used. Using the actual functional form of -Rtt.aa JFaurobertl (11987. 



1988) demo nstrated that Kneer's approximation used by iRees fc Salibal ( 119821 ) and ISaliba 



19851 . Il986l ) cannot account for the phenomenon of diffusion in frequency which takes place 
in the line wings. She showed that Kneer's approximation leads to non-negligible errors in the 
wing polarization obtained from slabs of finite thickness. In the case of a semi-infinite atmo- 
sphere this approximation affects both the line core and line wing polarization, as well as the 
intensity itself. Furthermore, she showed that in the presence of a background continuum, 
Kneer's approximation predicts incorrect values for both the position and the magnitude of 
the polarization maximum in the wings. IFaurobertl (ll988f ) considered also the influence of e 
and of the line strength on the linear polarization profiles and gave a qualitative interpreta- 
tion of the behavior of the emergent polarization profiles in terms of an Eddington-Barbier 
expression. 



The collisional redistribution ma trix given by iDomke fc Hubenvl ( If 9881 ) was used in po 



larized l ine transfe r computatio ns by iFaurob ert- Scholll ( 119921 . in plane-para llel atmospheres 



and b y iNagendral (119941 Il995l . in spherically-symmetric atmospheres). IFaurob ert- S choll 



(1992 



showed that th e line-wing polarization is sensitive to the elastic collisional rate Te 



while INagendral (119941 ) emphasize d that the l ine core po larization is sensitive to the depo- 
larizing collisions parameter . INagendral (119941 . Il995l ) presented a detailed study about 
the influence of e, the damping parameter a, the elastic collisional rate Te and on the 



Hubenv 


1985a: 


Frisch 


1980) 



In all the above-mentioned reference s, Feautrie r 's me t hod w as used to solve the PRD 
polarized line transfer equation, except in INagendral (11994 Il995[ ) who used a discrete space 
method. As is well-known, these type of numerical methods are computationally expen- 
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sive. Over the last two decades fast iterative methods based on operator-perturbation have 
been developed to solve the PRD radiative transfer problem of resonance line polarization 
ass uming a two-level model atom without lower- level atomic polarization (see the r eview 



by lNagendra & Sampoornall2009l ). More recently. ISampoorna fc Trujillo Buenol (120101 ) have 
developed very efficient and accurate symmetric Gauss-Seidel and Successive-Overrelaxation 
iterative methods to solve the above-mentioned PRD problem. Applying this fast radia- 
tive transfer method, in this paper we study in detail the effects of various atmospheric 
parameters on the linear polarization profiles formed under PRD conditions, but using semi- 
empirical models of the solar atmosphere and paying particular attention to understanding 
the shape of the emergent fractional linear polarization profiles in terms of the radiation field 
anisotropy within the medium. In particular, we consider two one-dimensional (ID) models 



of the quiet solar atmosphere: one based on the VALC model of IVernazza et al.l (1l98lL which 



has a relatively hot lower chromosphere) and the second based on the MCO model of lAvrett 



( 119951 . which has a relatively cool lower chromosphere). Although the solar atmosphere is 
highly inhomogeneous and dynamic and such ID models can only be considered as illustra- 
tive of the complex atmospheric conditions, they are suitable for achieving the main aim of 
this paper, namely a basic investigation on the impact of PRD effects on the shape of the 
emergent Q/I profiles and their sensitivity to various atmospheric parameters. In particular, 
we study the sensitivity of the linear polarization profiles to the following relevant quantities: 
(a) the photon destruction probability per scattering event e, (b) the strength of the spectral 
line, parameterized as r (see the discussion below Equation ([8]) for the definition of r), (c) 
the damping parameter a of the line absorption profile, and (d) the elastic collisional rate 

r E - 

The paper is organized in the following manner : in Section 2 we present a historical 
background on the PRD matrix for non-magnetic resonance scattering in a two-level atom. 
The formulation of the problem is presented in Section 3. A detailed study of the sensitivity 
of the linear polarization profiles to various atmospheric parameters and to the thermal 
structure of the solar atmosphere is presented in Section 4. Concluding remarks are given 
in Section 5. 



2. Redistribution Matrix for the Non-magnetic Resonance Scattering Problem 



In this section we present a brief backgro und on the PRD matrix for resonan ce line 
polarization. We follow the same notation as in ISampoorna &: Trujillo Buenol (l2010l ). 



The two-level atom PRD problem without polar i zation was formulated by iHummer 



(119621 ). Later, starting from the work of lOmont et al.l (Il972l ). more general PRD functions 
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includin g a better treatment of co 



rive d byj Heinzej (1198 



lision s , for both resonance an d subordinate lines, w e re de- 



fl.lHubenvl (Il982h . iHubenv et al.l f fl983al Jbh. Irlubenv fc Cooped Jl986h 



and iHubeny fc Lited (119951 ) . Revie ws on th e unpol arized PRD p roble m and its application 



to astrophysics were presented by IHubeny! (Il985bl ) and iFrischl (119881 ) . For a more recent 
review on the same topic , but w ith emphasis on numerical methods and radiative transfer 
modeling, see lUitenbroekl (120031 ) . 



When the polarization state of the spectral line radiation is taken into account the un- 
polarized PRD functions become 4x4 redistribution matrices that describe how the Stokes 
vector is redistributed in both frequency and angle. In the non-magnetic case the redistri- 
bution matrix, once averaged over the Maxwellian distribution of the scattering a toms, can 
be exp ressed as a product of the angle and frequency dependent PRD function of [Hummer 
( 119621 ) times a frequency-independent 4x4 phase matrix. Clearly, there is a intricate coupling 
between the frequency, angle and polarization state of the radiation field. 

A 4 x 4 phase matrix for t he J\ = — > J u = 1 — > Ji = scattering transition was 
derived by IChandrasekharl ( 119501 ) using classical ele ctrodynamics, w hich is referred to as the 
Rayleigh phase matrix. Using quantum mechanics iHamiltonl (119471 ) derived a more general 
Rayleigh phase matrix for a J\ — > J u — > Ji transition with arbitrary values for the lower 
and upper level angular momentum quantum numbers. For the particular case of a ID 
atmosphere, the 4x4 phase matrix reduces to a 2 x 2 phase matrix, due to the azimuthal 
symmetry of the problem. 

In polarized radiative transfer, using redistribution matrices that have intricate cou- 
plings bet ween frequency and a ngles is numerically expensive. Hence, to reduce the numer- 
ical work, iRees &: Salibal ( 119821 ) introduced the so-called hybrid approximation, in which it 
is assumed that the 4x4 redistribution matrix, depending on directions and frequencies, 
can be factorized into a 4 x 4 matrix, depending only on directions, times a scalar function 
depending only on frequencies (the so-called AA redistribution function). Indeed such a 
factorization is valid only in the atomic frame. In the laboratory frame it is mainly used 
as a practical "ansatz" for avoiding extremely time-consuming calculations and is justified 
only by means of heuristic arguments. Under the hybrid approximation, in a plane parallel 
atmosphere, the redistribution matrix can be written, in general, as a linear combination of 
terms of the form0 



1 Note that the only difference with Equations (33) and (37) of ISampoorna k, Truiillo Buenol (hold) is 

instead of g^ x ,, in order to be 



that in Equations ([1} and ([3]) of the present paper wg have written 4>xQ 
consistent with the definition of R given below. 



k 

xx' 
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R(x, x'; /i, //) = (1 - e) ^ P(/x, //) 



11 



Here x, x' are the non-dimensional frequencies of outgoing and incoming photons, and /i, ^ 
are the cosine of the polar angles 9, 8' with respect to the vertical direction in the atmosphere. 
R(x, x'\ /i, /i') gives the joint probability of absorbing a photon with frequency x' and direction 
fi' and re-emitting by spontaneous de-excitation a photon with frequency x and direction 
\i. The photon destruction probability is given by e = IV^! + T R ) with Y\ being the 
inelastic de-excitation collision rate and the radiative de-excitation rate. In Equation ([I]), 
9xx' = Rk.A A(x,x')/<i) x , whe re Rk,AA (with k — I, II, and III) are the AA redistribution 
functions of iHummerl (119621 ) and (p x is the normalized Voigt profile function. The function 
-Rfc,AA/0x gives the probability of absorption at frequency x', per emission at frequency x. 
Note that in Equation ([T]) we introduce for convenience the quantity g xx , which has the 
advantage of formally simplifying the expression for the line source vector (see Equations 
and ([TO]) below). Finally, for the case of a ID atmosphere, P(/x, //) denotes the 2x2 phase 
matrix. 

The redistribution matrix given in Equation (TjQ) does not take into account the elastic 
collisions, quantified by the parameters and D^ 2 \ However, the elastic collisional rate Te 
can be included into the hybrid approximation by replacing g xx , by a linear combination of 
type II and type III PRD functions, weighted by the factor 7 = (r^ + Ti)/(rR + Tj + Te), 
namely 

R(x, x'; fi, fi') = (1 - e) 0, [ 7 <& + (1 - 7)^] P(/i, //')■ (2) 

We remark that the (1,1) element of R given in Equations ([T]) and (J2J) is normalized to 
(1 — e), because it does not tak e into account the depolarizing collisional rate D^ 2 \ We 
refer the reader to lOmont et al.l (119721 ) for a detailed discussion on t he normalizat ion of 
the redistribution matrix. To qualitatively model the solar Ca II K line ISalibal (119851 ) used 
the redistribution matrix given by Equation but with a further simplification which 
consisted in replacing -Ri^aa by Kneer's (1975) approximation and Rm,AA by the CRD 
expression. Also, he neglected Ti in the expression for 7 (see above), and assumed that the 
matrix P(/x, //) wh ich mu l tiplies -Rhiaa (see Equation (j2J)) is the isotropic scattering phase 
matrix. Moreover, ISahbal (119851 ) did not take into account the depolarizing collisional rate 
£)( 2 ) i n his modeling. 



The redistribution m atrix given by lDomke fc Hubenyl (119881 ) for the unmagnetized case 
(see also lBommierlll997bl ) takes into account the effect of depolarizing collisions through the 

parameter. It is worth noting that this collisional redistribution matrix is yery g eneral, 
namely, it depends on the angle- dependent re distribution functions of lHummerl (jl962j). How- 
ever, for computational simplicity, following iRees fc Salibal (119821 ). some authors have used 
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the A A version of this gen eral collisional red istribution matrix (e.g., iFaurobert-SchollI 11992 
Nagendralll994l ) . Following iBommierl ( 1l997bl ) the A A version of the collisional redistribu tion 
matrix can be written as (see also Equation (37) of ISampoorna fc Trujillo Buenoll2010l ) 



R(x, x'; fi, //) = Yjk=o,2 W k{JuJu) 
where the branching ratios a and ft^ are given by 



a 



R 



r R + Ti + Ti 



(3) 



(4) 



ft 



R 



(5) 



r R + r : + £>(*) ' 

Note that = 0, and also that the factor (1 — e) is contained in the branching ratios. 
The coefficient Wo(Ji, J u ) — 1 ; and W2(Ji, J u ) characterizes the maximum linear polarization 
that can be produced in the line. In the case of a normal Zeeman triplet (J; = 0, J u — 1), 
W?.{Ju J„) = 1. Furthermore, P R (/x, ti') are the Rayleigh phase matrix multipolar compo- 
nents (ILandi DeglTnnocentil Il984l . see also Equation (35) of Sampoorna & Trujillo Bueno 
2010). 

The CRD limit is obtained from Equation ((3]) by taking Te 3> Tr (in other words a 1) 

and q ™, = 4> x > (compare the resulting expression with Equation (10.54) of lLandi Degrinnocenti fc Landolfi 
2004 . with the magnetic field strength set to zero). Notice that CRD means that there 
is no frequency correlation between the incoming and outgoing photons. In other words, 
the incident radiation is completely redistributed in frequency, so that there is complete 
non-coherence between the incident and scattered photon frequencies. This CRD limit is 
reached when the atomic system is i lluminated by a spectrally flat radiation field (e.g., 
Landi Degl'Innocenti &: Landolfi! 120041 ). Another good approximation to this CRD limit oc- 
curs when the atoms are so strongly perturbed by elastic collisions during the scattering 
process that t he excited elec trons are randomly redistributed over the substates of the up- 
per level (see iMihalasI Il978[ ). However, in the context of polarization, adopting the CRD 
approximation appears to lead to the following contradiction : CRD is obtained in the limit 
of very high collisional rates with respect to the radiative rate, but then the spectral line 
would be completely depolarized by depolarizing collisions. This apparent contradiction can 
be clarified by recalling the explicit physical meaning of the branching ratios appearing in 
Equation ([3]). First of all, for the problem at hand, elastic collisions produce two different 
effects represented by the collisional rates Te and D^ 2 \ The collisional rate Te is responsible 
for line broadening and for destruction of correlations between incoming and outgoing pho- 
ton frequencies. The rate is responsible for the destruction of the upper- level atomic 
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alignment and, thereby, causes depolarization. The branching ratio a gives the probability 
that a radiative decay from the excited state occurs before any type of collision (elastic or 
inelastic). Therefore, a gives the fraction of the scattering processes that are coherent in 
the atomic rest frame (i?n). The branching ratio (3^ gives the probability that radiative 
decay of the excited state occurs without destruction of the 2i^-multipole moment. There- 
fore, the total branching ratio \^ K ' — a\ gives the probability that radiative de-excitation 
occurs after a elastic collision (T E ) which redistributes the radiation in frequency without 
destroying the 2i^-multipole moment. Clearly \/3^ K ' — a\ gives the fraction of the scattering 
processes that are completely non-coherent in the atomic rest frame (Rm). Therefore, to 
get the CRD limit we only assume that elastic collisions that completely redistribute the 
radiation in frequency are sufficiently strong (Tg ^> Tr), so that the branching ratio a — > 0. 
In other words, and are assumed to be independent, even though their physical 
origin are the elastic collisions. In this case Equation (J3J) reduces to 

Rckd(z,s';a*,aO = Y,k=q,i W k{Ju Ju) 

x/^V^PfGuy), (6) 

where we have further replaced g 1 ^, by 4> x i. Clearly the branching ratio f3^ K > allows only that 
fraction of the scattering processes which do not lead to destruction of the 2ii"-multipole 
moment. 

We compute the -Rh,aa and -Rih,aa functions by numerically integratin g the correspond- 



ing angle- dependent redistribution functions (see Equations (59) a nd (Ql) oflBommierlll997bf) 



over all the scattering angles (see Equations (103) and (104) of iBommierl Il997bl ) . For the 
numerical integration over all the scattering angles we have used a 15 point Gauss-Legendre 
quadrature. In this paper we assume that the redistribution functions -Rh,aa and -Rhi,aa and 
the Voigt profile function <p x are constant throughout the atmosphere. Since these functions 
depend on the damping parameter a and on the reduced frequency x, which is defined as 
x = { u o — v)/Ai>y), with uq the line center frequency and Ah>n the Doppler width, the above 
assumption implies that a and Apd are also constant within the entire atmosphere. 



3. Formulation of the Problem 

In this section we discuss all the physical ingredients required for our investigation 
(presented in Section 4) about PRD scattering polarization profiles. Section 3.1 describes 
the basic equations of the problem. The model atmospheres used for our study in Section 4 
are discussed in Section 3.2. The definition of the anisotropy factor for PRD is considered in 
Section 3.3, where we also discuss several known facts about the anisotropy and its relation 
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to the source function gradient. In Section 3.4 we present the Eddington-Barbier relation 
for PRD problems, and show the direct dependence of the emergent polarization on the 
anisotropy within the model atmosphere under consideration. 



3.1. The Basic Equations 



We consider the standard two-level atom resonance line polarization problem in a one- 
dimensional, plane-parallel, static stellar atmosphere. Furthermore, we assume an unpolar- 
ized background continuum with no continuum scatt ering. The basic equations for t he ab ove- 
mentioned problem have been presented in detail in ISampoorna fc Trujillo Buenol (120101 . see 
their Section 5). Here we follow the same formulation, but for the sake of clarity we recall a 
few important basic equations from that paper. 



As shown by I Chandr asekharl (119501 ) , an azimuthally symmetric polarized radiation field 



can be described by the two Stokes parameters I XfM and Q x ^. The Stokes I Xfl parameter is the 
specific intensity, while Q Xfl quantifies the linear polarization (i.e., the difference between the 
intensity components parallel and perpendicular to a given reference direction in the plane 
perpendicular to the direction of the ray under consideration). In this paper, the positive Q Xfl 
direction is defined in the plane containing the direction of the ray and the vertical Z-axis. 
When polarization is taken into account, the Stokes source vector depends not only on the 
frequency x but also on \x = cos 8, with 9 the angle between the ray and the vertical Z-axis. 
It is possible to transform from the Stokes basis to an irreducible basis where the source 
vector components depend only on the fre quency x. Su ch a transformation is referred to as 
the "decomposition" of the Stokes vector. iFrischl (120071 ) has given a si mple way of achieving 
this g oal using the irreducible tensors for polari metry introduced by lLandi Degl'Innocenti 
( 1984 ). This decomposition technique was used in Sampoorna &: Trujillo Buenol ( 201oh . and 
we adopt it here also. 

In the reduced basis the irreducible intensity {I x ^q with K = and 2 satisfies the 
following transfer equation 

d 



K i 



dr. 



(7) 



Xft, 



where the irreducible source vector components are given by 

sK MSix)o+rBU? 



OS 



xJO 



_|_ r 



Here r Xil is the total optical depth defined by o\t xim = — {xi^x + Xc)dz/ fi, with z the distance 
along the normal to the atmosphere, xi an d Xc are the line and continuum opacities, r = 
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Xc/Xh B i s Planck function, and Uq is equal to unity for K = and zero for K = 2. 
The irreducible line source vector components are given by 



(S lx )£ = eBU« + W K {J h J„)(Jx)o . (9) 

where 

/+oo 
(io) 
-oo 

The angle integrated irreducible tensors of the radiation field are respectively given by 



1 



-i 



« = 2 J dfxl xfl , (11) 

{J*)l = ^= £ d/i [(3/i 2 - 1)1^ + 3(/i 2 - 1)Q,J . (12) 

The irreducible components of the radiation field and of the source function are connected to 
the Stokes p arameters and to the source vector through a simple formula (see Appendix B of 



3_Dc 

Frischl 120071 ). Note that in the solar atmosphere the dominant contribution to (J x )l is given 



by the first term in the square bracket of Equation ( |T2l) . 



3.2. The Model Atmospheres 

In this section we describe the two solar model atmospheres that we have chosen for our 
detailed study of the linear polarization profiles presented in Section 4. As already mentioned 
in Section 1, we consider the VALC and the MCO models for this purpose. Figure Q] shows 
the variation of the temperature as a function of height for the VALC (dotted line) and 
MCO (solid line) models. We note that both models include the transition region as well as 
the corona. Most of the results presented in this paper correspond to spectral lines which 
become optically thin for heights z > 2100 km. However, we also present some results for a 
very strong spectral line whose line center originates in the transition region of the VALC 
model. 

For a hypothetical line at A = 5000 A we calculate the Planck function B (in ergs -1 cm -2 sr -1 cm -1 ) 

afl 

2hc 2 1 

B = e hc/(Xk B T) _ l ' ( 13 ) 



2 Rigorously speaking, the Planck function expression should be that corresponding to the Wien limit 
because we are neglecting stimulated emission. However, for visible and UV lines both expressions give 
similar values. 
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where c is the speed of light, h is the Planck constant, and fee is the Boltzmann constant. As 
the tabulated height grid of the VALC and MCO models is crude, we actually first interpolate 
the tabulated temperatures on a much finer height grid. We use cubic spline interpolation for 
this purpose. Our height grid has a spacing of Az = 5 km. Note that achieving the accuracy 
corresponding to such a fine spatial grid requires finding the self-consistent solution via the 
application of a highly- convergent it erative scheme, such as that on whic h the PRD computer 



program used here is based on (see ISampoorna fc Trujillo Buenoll20!0l ). 



We assume that both solar model atmospheres are exponentially stratified, so that the 
continuum optical depth is described by the law r c = 2.2 exp(— z/H) with z (in km) the 
vertical height in the atmosphere, and H the scale height (which we have chosen equal to 
120 km). Moreover, we assume that the line strength parameter r, the photon destruction 
probability per scattering e, and the damping parameter a of the line absorption profile are 
constant with height. Figure |2] shows a plot of r c versus height for the V ALC model. For 



r c we use the tabulated value at A = 5000 A from IVernazza et al.l (119811 ) . Overplotted is 
our continuum optical depth versus height. Clearly, the r c of the VALC model does not 
show a single scale height. In fact, our attempts to fit the r c of VALC required the use of a 
weighted combination of two exponentials with two different scale heights. We find that for 
heights lower than 600 km one has to use a scale height of 60 km, while for larger heights one 
needs a scale height of 250 km, in order to reasonably fit the r c of VALC. However, to keep 
the problem as simple as possible, we choose a single scale height of 120 km for the entire 
atmosphere and Tq = 2.2 that nearly coincides with the r c value of the VALC model at the 
lower boundary of the atmosphere (see dotted line in Figure [2]). The same choice applies 
to the MCO model. This strategy appears to be reasonable for facilitating a systematic 
investigation on the impact of PRD effects on the shape of the emergent Q/I profiles. 



3.3. The Anisotropy 

The (J x )o tensor, which is dominated by the contributio n from the Stokes I Xjl parameter, 



chara cterizes the 'degree of anisotropy' of the radiation field ( jLandi Degl'Innocenti fc Landolfi 



2004 . see also Trujillo Bueno 2001). We define a 'monochromatic anisotropy factor' as 

A = {J x )l/{J x )l (14) 

For clarity of notation, we denote the monochromatic anisotropy A for CRD and PRD by 
Acrd and Aprd, respectively. We also define a 'frequency integrated or mean anisotropy 
factor' for PRD as 

^PRD = (Jx)o/(Jx)o, (15) 
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where (J x )o and {J x )q are defined as in Equation ffTOj) . In the case of CRD the frequency 
integrated anisotropy factor is defined by 



Clearly, the mean anisotropy Acrd does not depend explicitly on frequency. However, in the 
solar atmosphere Acrd changes with height. It then results that, even in CRD, since different 
points of the profile form at different heights, we have an "apparent" frequency dependence 
of Acrd across the profile. On the contrary, Arrd depends explicitly on frequency and this 
effect adds to the previous one. 

In Figure [31 we illustrate the important difference between the monochromatic and mean 
anisotropy factors for both, the PRD and CRD cases. To compute these anisotropy factors 
we have used the VALC model atmosphere (see Section 3.2) and e = 10~ 4 , r = 10 -5 , a = 10~ 3 
and r E /r R = /r R = 0. In Figure [3] we plot A and A versus the reduced frequency at 
the height where for a line-of-sight (LOS) with \x = 0.11 the condition r Xfl = 1 is satisfied. 
Clearly, the monochromatic anisotropy factors Aqrd and Aprd are similar. The differences 
are basically due to the slight differences that are found in Stokes I Xfl for CRD and PRD 
(compare the dotted lines in the intensity panels of Figures H] and [5]) . It is interesting to note 
that Ackd(t x ^ = 1) does not tend to zero in the wings. This is because the condition r XfJi = 1 
for a LOS with \x = 0.11 is satisfied at the height of 360 km for all frequencies x > 20. At this 
height the anisotropy Acrd for x > 20 is determined by the Planck function gradient at that 
height. The mean anisotropy factors Aprd and Acrd differ greatly, particularly in the wings. 
This can be understood as follows. In the case of CRD we integrate {J x )1 and (J x )q over all 
x after multiplying them by the absorption profile function <p x (see Equation (j!7p ). which 
goes to zero in the wings. Thus, at any depth T XfM the dominant contribution to ( Jq )crd 
and, thereby to Ackd comes from the line core region. Furthermore, the monochromatic 
anisotropy Ackd decreases with depth, particularly in the line core. As already noted above, 
the frequency dependence of A C kd is due to the fact that different parts of the line profile 
are formed at different heights. Combining these facts, it is easy to see that A C kd at r X(U = 1, 
decreases as the frequency increases. In particular, far in the wings of strong lines like the one 
we are considering in Figured (i.e., with r = 10 -5 ) Ackd is zero (even if Ackd is substantial), 
because we are looking deep in the atmosphere where Ackd in the line core is zero (which 
dominantly contributes to A CRD ). In the case of PRD, we integrate (J x )q and (J x )q over all 
x after multiplying them by the type II redistribution function (see Equation fllOp ). which 
behaves like CRD in the core and like coherent scattering in the wings. Thus, unlike CRD, 



^CRD — (^o)cRd/(^o)cRD, 



(16) 



where 




(17) 
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at any depth r Xfl the mean anisotropy Aprd depends explicitly on x. Furthermore, in the 
line core it looks somewhat similar to A CRD . But in the wings A PRD coincides with ^4p R o, 
due to the coherent behavior of g^ x ,. 

The importance of the anisotropy factor, which is the fund amental quant i ty tha t de- 



termines the shape of the emergent polarization, was realized by iRees fc Salibal ( 119821 ) who 
pointed out that the polarization is a mapping of the depth dependence of the anisotropy of 
the radiation field within the atmosphere. However, they neither quantified nor presented a 
detailed stud y of the anisotropy factor. Some information can be found in the PhD thesis of 



Salibal (ll986f ). who presented a study of the anisotropy for the CRD case (see his Chapter 4). 
He also illustrated the relation of the anisotropy factor to the source function gradient by 
considering, in a finite slab atmosphere, a simple piece-wise linear source function. His quan- 
tities a and a are quite similar to our quantities A CRD and A CRD . However, in the case of 



PRD, ISalibal ( 119861 ) did not define a mean anisotropy factor Ap R D similar to that defined in 
our Equations (|T5|) and ( TTOj) as he used Kneer's approximation for -Ri^aa- Nevertheless, he 
correctly noted that for line core frequencies the anisotropy of the radiation field is controlled 
by ^4crd and for the line wings it is controlled by Appj> In this paper we show that a better 
way of defining the anisotropy for the PRD problem is through Equation (I15p along with 
Equations (fT0|)-(fT2l). 

A Milne-Eddington atmosphere is a very suitable model to clearly demonstrate the 
strong dependence of the anisotropy fa ctor on the gradient of the source function for the 



Stokes I x „ parameter (see Section 3 of iTrujillo Buenol 120011 ). The mean anisotropy A C rd 



is bounded by —1/2 < \/2^4 CRD < 1, the lower and upper bounds being for a purely hori- 
zontal radiation field without azimuthal dependence and for a purely vertical illumination, 
respectively. We summarize the following facts about the anisotropy in a stellar atmosphere : 

(1) The anisotropy is essentially negative for an atmosphere with no or very small gradient 
in the source function. The negative values of the anisotropy means that the radiation field 
is predominantly limb-brightened (i.e., predominantly horizontal). 

(2) The larger the source function gradient the larger the anisotropy. Positive anisotropy 
values imply that the radiation field at the spatial point under consideration is predominantly 
limb-darkened (i.e., predominantly vertical). 

In summary, in a stellar atmosphere the outgoing radiation is predominantly vertical 
(which makes a positive contribution to ( J x )o), while the incoming radiation is predominantly 
horizontal (which makes a negative contribution to (J x )q). More precisely, "vertical" rays 
(i.e., with \jx\ > l/v^3) make positive contributions to A, while "horizontal" rays (i.e., with 



u| < 1 / \/3) make negative contributions to A. As seen in Figure 4 of ITrujillo Bueno 



(120011 ). the larger the gradient of the source function the greater the anisotropy factor of 
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the pumping radiation field, and the larger the amount of atomic level polarization. A more 
detailed stu dy of the anisotropy factor and i t s rela tion to the source function gradient can 
be found in lLandi DegFInnocenti &: Landolnl (120041 ) . 



More recently, iHolzreuter et al.l (120051 ) have presented a study of the monochromatic 
anisotropy Aprd for the solar Ca I 4227 A and Nal 589 A lines taking into account the AA 



collisional redistribution matrix (see, e.g., lStenflolll994j ). They also illustrate the relation of 



^4prd to the source function gradient. However, as they consider only A PRD the interpretation 
of the emergent linear polarization in terms of Appj) is not straightforward. In this paper we 
show that a better quantity that can be directly related to the emergent linear polarization 
is v4 PRD . 



3.4. The Eddington-Barbier Relation 



In the case of CRD, the emergent polarization directly depends on the anisotropy of 
the radiation field (see Equation (13) of iTrujillo Buenolll999l . see also Section 1 of Trujillo 
Bueno 2003). Such a relation can be considered as the generalization of the Eddington- 
Barbier relation for scattering polarization. Thus, the emergent polarization is governed by 
the anisotropy of the radiation field, which in turn is strongly related to the gradient of the 
source function for Stokes I Xfl and, t herefore, to th e jx dependen ce of I x ^. Note that the 



Eddington-Barbier relation given by ITrujillo Buenol ( 119991 . 120031 ) for the CRD problem is 
very general, as it takes into account the effects of lower-level polarization. 



An Eddington-Barbier relation for the PRD case was given by iFaurobertl ( 119871 . Il988[ ) 
who used it for a qualitative analysis of the linear polarization profiles formed in an isothermal 
model atmosphere. According to the Eddington-Barbier relation the emergent intensity 
and the emergen t polarization from a semi-infinite atmosphere can be written as (see also 
Equation (41) of lFrisch et aDl2009h 



0) ~ (S x f (r x , = 1), 
Qx,(r xll = 0) ~ ^Gu 2 - l)(S x ) 2 (r xp = 1) 



18) 



Restricting to the pure line case, and neglecting the thermal term eB, it is easy to show from 
Equations (19"])- (1T2"]) and (USD, that 



Q 



Ixh 



' X/l 



0) 



2V2 



(fi 2 -1)W 2 A prd (t c 



1). 



(19) 



This approximate formula shows clearly that the emergent fractional linear polarization 
depends on the mean anisotropy Aprd at t x/1 = 1 along the LOS. In this paper, we therefore 
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show Aprd versus the reduced frequency at the height in the model atmosphere where r Xfl = 1 
for a LOS with /i = 0.11. Furthermore, as our reference direction for positive Stokes Q is 
defined perpendicular to the limb, we plot —Q/I whose positive values indicate polarization 
parallel to the limb. 

4. The Emergent Linear Polarization Profiles 

In this section we present a detailed study of the PRD scattering polarization profiles 
computed in the VALC and MCO model atmospheres (see Section 3.2 for details on the 
assumptions made). For these two atmospheric models we present our results by giving, for 
each choice of atmospheric parameters, a set of figures with four panels showing: (a) the 
emergent intensity at fi = 0.11 (normalized to the continuum intensity) as a function of the 
reduced frequency; (b) the emergent fractional linear polarization at fi = 0.11 as a function 
of the reduced frequency; (c) the mean anisotropy factor, A, evaluated at the height where 
the monochromatic optical depth along a LOS with fi = 0.11 is equal to unity, as a function 
of reduced frequency; (d) the line source function at line center (x = 0) for the intensity 
propagating along a LOS with fi = 0.11 (denoted as Si) as a function of height. For the 
sake of comparison, the corresponding CRD linear polarization profiles are shown for certain 
cases. 



4.1. Results of calculations in the VALC Model 

In this section we study the influence of the line strength parameter r, the photon 
destruction probability per scattering e, the damping parameter a, the elastic collisional rate 
T E and the depolarizing rate on the PRD linear polarization profiles. We use A = 5000 A, 
e = 10~ 4 , r = 10~ 5 , a = 10 -3 , and T^/Tr = /Tr = as the nominal case, around which 
we vary the various parameters mentioned above. Such a study should help us to understand 
the results of future, more realistic computations where all the above-mentioned parameters 
are depth dependent. Unless stated otherwise we set W 2 — 1, which represents the case of a 
line transition with Ji = and J u — 1. For the sake of clarity we divide this section about the 
results of our radiative transfer calculations in the VALC model in various subsections dealing 
with the sensitivity of the solution to r (Section 4.1.1), e (Section 4.1.2), a (Section 4.1.3), 
r E /r R and £> (2) /r R (Section 4.1.4). Some additional studies are presented in Section 4.1.5 
concerning lines with wavelengths other than A = 5000 A, lines formed in the transition 
region, and lines with other W 2 values. 
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4-1.1. Effect of the Line Strength Parameter r 

Figure H shows the response of I(r Xfl = 0), -(Q/I)(T Xft = 0), Aprd^ = 1), and 
Si(x = 0) for a LOS with \i = 0.11 to variations in r. For comparison the corresponding 
CRD case is shown in Figure [5j For notational simplicity we drop the subscripts xfi on / 
and Q. Table [T] gives the height in the VALC model atmosphere at which lines of different 
strengths have r 0A1 = 1, both for /i = 0.11 and /x = 1. 

We first discuss the PRD line profiles presented in Figure H] and then compare them with 
the CRD profiles of Figure [51 Table [1] indicates that as r decreases we progressively sample 
lines that are formed higher in the atmosphere. As r increases the line becomes weaker 
and the broad damping wings disappear (see Figure H]). The —Q/I profile of the line with 
r = 10 -6 shows a maximum in both the line core and in the line wing (i.e., at about « 15 
Doppler widths) as well as a minimum in the near wing (i.e., at about w 4 Doppler widths). 
As r increases both the wing maximum and the wing minimum are shifted towards line center 
and decrease in magnitude, to eventually disappear (see the long-dashed line in Figure HJ). 
The line core polarization, however, initially decreases and becomes negative for r = 10~ 2 
and then becomes positive again for r = 0.7. This response of —Q/I to the variations in r can 
be understood using the mean anisotropy Aprd plotted in Figure HI It is worth to note that 
the shapes of the —Q/I profiles that we obtain with PRD fo r r values between 10~ 6 and 10~ 2 



(see F igure HI) are similar to those classified as M-signals by lBelluzzi fc Landi Degl'Innocenti 



(120091 ). whereas the shape of the —Q/I profile for r = 0.7 is similar to that cla ssified as S- 



signal by these authors. We recall from iBelluzzi fc Landi DeglTnnocentil (120091 ) that a Q/I 



profile with a single peak at line center is classified as S-signal, while Q/I profiles showing a 
polarization maximum in the wings, a decrease in amplitude approaching the line core, and 
eventually a narrow peak at line center are classified as M-signals. 

For r = 10 -6 and 10~ 5 the —Q/I profile is nearly proportional to Appx> in the frequency 
range < x < 5 Doppler widths, so that Equation ( fl~9l) is applicable. However for x > 5, 
one cannot neglect the contribution of r in (S x )q as was done to deduce Equation (fT9l . In 
the presence of a background continuum, one can approximately write the emergent Q/I 
profile as (still neglecting the thermal term eB) 

T (T - " 0) ~ 271 (M " } 2 <t> x + rB/{J x )l ■ m 

The term with r, in the denominator of Equation ( 120|) . is responsible for the formation of a 
wing maximum in —Q/I. The fact that the introduction of r gives rise t o a w ing maximum 



in —Q/I was already noted for isothermal models by iRees fc Salibal (Il982[ ). This wing 



maximum is located at the frequency where the radiation field starts being influenced by 
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the continuous absorption (see lFaurobertlll988l ). In other words, at that frequency the line 



source function and the (unpolarized) continuum source function equally contribute to the 
total source function. For isothermal model atmosphere the p osition of this wing maximum 



can be estimated using the thermalization frequency given by iFrischl (119801 ) 



The frequency range where Equation (|T9|) can be applied decreases as r increases and can 
no longer be applied for r > 10 -2 . On the other hand, when r increases, the term containing 
r in Equation (1201) starts influencing not only the wings but also the line core region, thereby 
greatly reducing the polarization in the wings to zero and confining the —Q/I profile to the 
line core. It is worth to note that for r < 10~ 3 the mean anisotropy Aprd is positive in 
the line core and in the line wings showing that the radiation field is limb-darkened. For 
those frequencies x whose r Xfl = 1 at the heights where we have the VALC temperature rise 
(~ 500 - 1000 km; see Figured]), the anisotropy is negative, showing that the radiation field 
is limb-brightene d. This fact was already noted for schematic chromospheric-like models by 



Faurobertl (119881 ). 



As r increases, the anisotropy Aprd in the line core decreases and even becomes negative 
for r = 10 -2 . This is because the line core photons now start to originate in deeper layers of 
the atmosphere, where the radiation field for the line core photons tends to become more and 
more isotropic. For r = 10 -2 a limb-brightened radiation field dominates to give a negative 
Aprd in the line core. However, for r = 0.7, the anisotropy Aprd reaches a large positive 
value. Note that the corresponding line source function Sj shows a larger departure from 
the Planck function, not only above the temperature minimum, but also in the photospheric 
range of heights 300 - 500 km (where actually the line core is formed; see Table [1]), differently 
from what happens with other values of r (see the Si panel of Figure H]). In fact, by taking 
values of r between 10~ 2 and 0.7, it can be shown that there is a gradual transition in the 
line-center anisotropy A PRD from negative values to a large positive value. 

It is worth to note that Aprd in the wings attains a constant value, which is the same 
for all the values of r. However, the frequency distance from line center at which this 
constant value is reached decreases as r increases. This can be understood by taking the 
example of r = 10~ 5 . In this case we find that the Stokes parameters at all frequencies 
x > 20 are "formed" at a height of 360 km (in the sense that this is the height at which 
the monochromatic optical depth is unity). When we plot the total source function for all 
these frequencies we find that they nearly coincide with the Planck function and, moreover, 
that they are dominated by the continuum source function. Thus, the anisotropy at those 
frequencies is determined by the Planck function gradient at a height of 360 km. As r 
increases the frequency value above which all the photons form at 360 km progressively 
moves towards line center. 
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A comparison of Figures H] and [51 shows that CRD is a good approximation in the line 



core and parti c ularly at line center. This is a well known result sin ce the deta i 



Dumont et al.l (119771 using the type I redistribution function) and iFaurobertl (119871 . [1988 



ed work by 



using the type II redistribution function). Furthermore, for sufficiently large values of r 
(weak resona nce lines) the CRD approximation can be safely used to model line profiles 
( Frisch 1996), like for insta nce the Q/I profiles of the photospheric line of Sri at 4607 A 



sec 



Faurobert-Scholll 119931 ) . This is because these lines do not have well-developed broad 



wings (unlike the strong resonance lines). Therefore transfer of photons i n the line wings, 
wher e frequency coherent scattering can play a crucial role, is negligible (see iFaurobert-Scholl 
19931 ). As a result, even in the PRD case, the polarization gets confined to the line core as 



in the CRD case. 



4-1-2. Effect of Photon Destruction Probability Per Scattering e 

Figure E] shows the sensitivity of /, —Q/I, ^prd^^ = -0 anc ^ = 0) at // = 0.11 to 
changes in the e value. The corresponding CRD case is shown in Figure [71 

When e increases the coupling between the radiation field and the Planck function 
becomes stronger. From Figure El we clearly see that Si ~ B in a greater region of the 
atmosphere for larger values of e. As a result, the intensity increases. However, the anisotropy 
decreases in the line core. Moreover, for frequencies x ~ 2.5 the absolute value of the 
anisotropy (which is here negative) increases, while for x > 7 the anisotropy is independent 
of e. Note that x ~ 7 corresponds to the wing maximum in the —Q/I profile. As already 
noted in Section 4.1.1, at this frequency the total source function starts being dominated by 
the unpolarized continuum, so that the anisotropy as well as —Q/I are independent of e. A 
comparison with the corresponding CRD case in Figure [7J shows that the above discussion 
for the line core region can be used to understand the CRD profiles also. 



4-1.3. Effect of the Damping Parameter a 

Figure M shows the response of I, —Q/I, -Aprd = 1) and Si(x = 0) at \i = 0.11 to 
variations of a. The corresponding CRD case is shown in Figure [9j 

The larger the damping parameter a the broader the absorption profile. As a result the 
emergent / as well as —Q/I profiles broaden as a increases. Furthermore, an increase of a 
implies that the source function departs from the Planck function to a slightly larger extent 
(see Figures [8] and [9]) . In the CRD case the anisotropy decreases with increasing values of a 
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and the —Q/I profile does the same. In the PRD case the anisotropy at line center does not 
show a large sensitivity to a, and thereby the same happens with —Q/I. On the contrary, 
for x > 3, Aprd shows a large sensitivity to variations in a. In particular, the anisotropy 
increases in absolute value for frequencies that are formed in the temperature rise region of 
the VALC atmosphere. Furthermore, as a increases the frequency at which Aprd reaches 
a constant value also increases. Consequently the negative dip (or wing minimum) and the 
wing maximum in —Q/I are very sensitive to a. The position of the negative dip as well 
as that of the wing maximum shifts away from line center. Since the position of the wing 
maximum is determined by the frequency at which the total source function starts being 
dominated by the unpolarized continuum, it is clear that an increase in a results in a shift 
towar ds larger frequencies of the wing maximum of —Q/I. It is worth to note that ISaliba 
( 119851 ) also shows the influence of a on / and —Q/I, though a detailed analysis is missing. 
Nevertheless, it is clear from our Figure Elthat a height dependenc e of a h a s imp ortant effects 
on the emergent polarization profiles, a point already stressed by ISalibal (119851 ) . 

We point out that our assumption of a constant damping parameter a throughout the 
atmosphere is actually in contradiction with the hypothesis of an exponentially stratified at- 
mosphere. This is because, as a is proportional to the density, it should also be exponentially 
stratified. We also point out that the damping parameter a is proportional to the elastic 
collisional rate Te, that has been here assumed to be zero. Notwithstanding these inconsis- 
tencies, we considered worthwhile to perform our studies as if the damping parameter were 
an independent quantity. 



4-1-4- Effect of the Elastic Collisions 



Figure [10] shows the sensitivity of /, —Q/I, v4prd(t^ = 1) and S\(x = 0) at fi 
to the variations of Te/Fr. The corresponding CRD case is shown in Fi gure 



iflcj 



0.11 
The 



depolarizing elastic collision parameter is assumed to be = 0.5 T E (e.g., IStenflol Il994j ) . 
We recall that in the CRD case, the redistribution matrix is given by Equation Clearly 
only is the relevant quantity for CRD. Therefore, for CRD varying Te/Tr is equivalent 
to varying D^ 2 > /Tr since we have assumed D^/Tr = 0.5 Te/Fr. We note t hat Te/Fr mainly 
operates in the line w ing, while /Tr operates in the line core (see also iFaurobert-Scholl 
19921 ; iNagendral 1 1994 ) . In the solar atmosphere the elastic collisions Te and increase 
with depth in the atmosphere, while they are negligible in the upper layers, namely in the 
upper chromosphere and transition region. Thus, depth dependent elastic collisions affect 
only the wings through Te/Fr, while the line core i s nearly unaffected as it is formed higher 
in the atmosphere where is negligible (see also lFaurobert-Scholllll992l ). 
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As expected, elastic collisions have little effect on the line source function Si. In the 
case of PRD, as Te/Tr operates in the line wings, there is here some sensitivity to this 
parameter in the intensity profile. However, for CRD, as /Tr operates mainly in the line 
core, there is hardly any effect on the intensity. As Fe/Fr or /Tr increases, basically 
the mean anisotropy factors Acrd and Aprd decrease at all frequencies. In the case of PRD, 
with an increase in elastic collisions, the contribution of g™, increases. Since g™, behaves 
more like CRD, when its contribution increases Aprd i n the wings decreases and it then 
gradually tends to zero. 



4-1.5. Additional Studies 

Here we present some additional studies related to the PRD linear polarization profiles. 
From now on we do not present the corresponding CRD cases. 

In Figure [T2] we show J, —Q/I, Aprd(t x/< = 1) and Si(x — 0) at fi — 0.11 for two 
different spectral lines. Clearly, the gradient of the Planck function as well as of Si are larger 
at shorter wavelengths (e.g., for A = 3933 A compared to those for A = 5000 A), particularly 
in the photosphere. Therefore, when we move to the blue part of the spectrum the anisotropy 
increases particularly in the wings, so that the polarization also increases there. 

Figure [TBI shows /, —Q/I, Apkd(t xi1 = 1) and the line source function Si(x = 0) at 
/i = 0.11 for lines formed in the transition region of the VALC model. As r decreases, the 
line is formed higher in the atmosphere (see Tabled]) where it encounters a steep temperature 
rise, so that the radiation field tends towards becoming limb-brightened. As a result, the 
anisotropy decreases in the line core, while it increases in absolute value at the negative dip 
(3 < x < 7). Moreover, as r decreases Equation ffT9l) can now be applied to a larger frequency 
domain and hence the wing maximum in —Q/I, where the unpolarized continuum source 
function starts to dominate, now shifts towards larger frequencies. As the spectral line is 
formed in the region of the steep temperature rise, the intensity shows a typical emission 
profile, increasingly self-absorbed. 

Figure [U shows the effect of W 2 on -Q/I and A(r Xfl = 1) at // = 0.11, for the VALC 
model including its transition region and corona. Since W 2 has only marginal effects on the 
intensity and the line source function Si, we do not show them here. As expected, when 
Wi decreases the polarization decreases everywhere. However, it is interesting to note that, 
in general, we expect a linear dependence on W 2 (since (Si x )q directly depends on W 2 , see 
Equation (jUJ)). In other words, as W 2 is reduced to half, we expect the polarization to be also 
reduced to half its value throughout the profile when compared to the corresponding W 2 = 1 
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case. Instead, we see that for x < 3, —Q/I(W2 = 0.5) is nearly 0.5 x —Q/I(W2 = 1). But 
for x > 3 there is no linear dependence on W 2 . This can also be seen in the mean anisotropy 
^4prd, which does not contain the W 2 factor in its definition (see Equations (JT0]) - (fT2"]) and 
f fT5"j) ). If there were a linear dependence on W 2 , the mean anisotropy A PRD should have been 
exactly the same for both W 2 = 0.5 and W 2 = 1. Instead, we see considerable difference 
for x > 3. This difference can be attributed to the contribution of Q to (J x )o, as I is only 
marginally affected by the value of W 2 . Thus, we conclude that the scattering polarization 
de pends non- linearly on l^. This is true even for the CRD case. This fact was also noted 
bv lstenflo fc Stenholml Jl976h . 



4.2. Results of Calculations in the MCO Model 

Here we study the effect of a relatively cool lower chromosphere on the emergent Q/I. 
As the sensitivity to parameters like r, e and a is already quite well understood for the 
VALC model, here we only present : (1) the comparison between the line profiles of VALC 
and MCO models, and (2) the influence of the elastic collisional rates Ee/Tr and /T-&. 
We do not show the corresponding CRD the main effects are confined to the line 

core. 

Figure [TS] shows the effect of the thermal structure of the atmosphere on the polarized 
line formation. From Figure [1] we see that the temperature minimum is shifted from the 
height z = 500 km in the VALC model to 1000 km in the MCO model. Thus, monochromatic 
radiation formed in the region 500 - 1000 km in the VALC model is limb-brightened, while it 
is limb-darkened in the MCO model (as there is no temperature rise in that region). As the 
region of temperature rise for the MCO model is located between 1000 and 2000 km, where 
the line core is formed, we find that the anisotropy in the line core is smaller for the MCO 
model than for the VALC model. The temperature rise tends to make the radiation field 
more limb-brightened. Finally, in the wings the anisotropy is smaller compared to that of 
the VALC model. This is because frequencies x > 4 are formed in the region 360 - 500 km, 
where the temperature gradient is larger for VALC than for MCO. 

Figure [16] shows the effect of varying the elastic collisional rate Fe/Tr for the MCO 
model. Effects of varying r E /r R or D (2) /r R for MCO model is more or less similar to that 
for the VALC model. We note that when T^/Tr increases, the differences in the anisotropy as 
well as in the polarization obtained from the MCO and VALC models decrease (for example 
compare the long-dashed lines in Figures [TOl and [T6|) . 
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Conclusions 



The interpretation of the second solar spectrum is an exciting and tough challenge in 
Solar Physics. To achieve this it is essential to solve the polarized radiative transfer equation 
taking into account various physical phenomena. In the case of strong resonance lines one of 
the important physical ingredients required for the modeling of the linear polarization profiles 
is the so-called partial frequency redistribution (PRD). In this paper, therefore, we have 
studied the scattering polarization profiles forme d under PRD. We use the angle-average d 
version of the collisional redistribution matrix (see lDomke fc Hubenylll988l ; lBommierlll997bl ). 
which is valid for a two-level atom with unpolarized lower level. For the numerical method 
of solution we have applied the very efficient and accurate symmetric Gauss-Seidel method 
presented in ISampoorna &: Trujillo Buenol (120101 ) . 



Our main focus in this paper has been to study the influence of various atmospheric pa- 
rameters (namely, e, r, a, Te/Pr, D^/Tr and the temperature structure of the atmosphere) 
on the linear polarization profiles. Interestingly, for strong lines the —Q/I profiles that we 
obtain with PRD, for a wide range of the various parameters, are o f the shape that has 
been classified as M-signals by iBelluzzi fc Landi DeglTnnocentil ( 120091 ). We recall that M- 
signals represent Q/I profiles with a three-peak structure. Therefore, we find that for strong 
resonance lines the three-peak structure for Q/I is a common feature for various combina- 
tions of the other relevant parameters, thereby confirming previous results either obtained 
using more crude model atmospheres (isothermal atmospheres) or with less sophisticated 
modeling. 

We interpret the linear polarization profiles using the Eddington-Barbier relation and 
the anisotropy at the atmospheric height where r Xfi = 1. To this end, as in the CRD 
case, we define a suitably frequency integrated or mean anisotropy factor A for PRD, that 
can be directly related to the emergent linear polarization. We have shown in great detail 
that the emergent scattering polarization is a mapping of the anisotropy at r Tjl , = 1 (cf., 



Rees & Saliballl98a iTruiillo Buenoll200ll : lLandi Deri'Innocenti fc Landolfill2004f l. Our study 
shows that a height dependent damping parameter a has significant influence on the emergent 
linear polarization profiles, particularly in the wings of strong resonance lines. In general, 
scattering polarization profiles are very sensitive to the thermal structure of the atmosphere 
(see Figure ITS]) . For the range of parameters considered in this paper, we find also that the 
dips in Q/I are located at values of reduced frequency x ranging from 2 to 5. Furthermore, 
the values of Q/I at these dips (including their sign) seem to be very sensitive to the thermal 
structure and to the various atmospheric parameters. 

The results of our investigation emphasize that PRD is required to model the linear 
polarization of strong resonance lines like Ca II K and Mg II k (compare Figures H] and EJ see 
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also Figure 113]) , while CRD is a good approximatio n to model weak resonance lines like Sri 
4607 A (see also lFaurobert-Scholllll993HFrischlll996l . and compare the long-dashed lines of our 
Figures S] and . A strong resonance line is characterized by broad dam ping wings, where 



frequency coherent scattering can play a crucial role (e.g., iFrischl Il996l ). Weak resonance 
lines lack such damping wings and, therefore, the corresponding linear polarization, even 
computed in PRD, gets confined to the line core, with a single peak at line center, similarly 
to what is obtained when assuming CRD. In order to quantify the difference between the 
linear polarization profiles computed with CRD and PRD for the case of weak resonance 
lines, we define the following two quantities : 



5 



Q 



(x = 0) 



|(Q//)prd(x = 0)-(Q/I) C rd(x 
|(Q//)prd(x = 0)| 



(21) 



which gives the relative difference in amplitude of Q/I computed using CRD and PRD at 
line center, and 

|(HWHM) PRD - (HWHM)crd| 



5 [HWHM] 



(22) 



|(HWHM) PRD | 

which gives the relative difference in the half width at half maximum (HWHM) of the Q/I 
profiles computed using CRD and PRD. Taking the particular case of r = 0.7 (see long- 
dashed lines in Figures H] and [5]) , we varied e, a, and Te/Pr as done in Figures 151 - [TT1 for 
r = 10" 5 . We find that 5 [(Q/I)(x = 0)] is around 5% and 5[HWHM] is about 10% for the 
above-mentioned parametric study with fixed r = 0.7. Clearly, the error introduced through 
the use of CRD for weak resonance lines is quite small. Therefore, we may conclude that 
CRD is a good approximation for weak resonance lines. 

Future investigations should extend the present work by considering angle-dependent re- 
distribution functions and the Hanle effect of deterministic magnetic fields. On the one hand, 
radiative transfer cal culations with angle-dependent redistri bution functions are numerically 
very expensive (e.g.. lFaurobertlll988l ; iNagendra et al.ll2002l ). To assume the angle-averaged 



PRD function used here instead of the general a ngle-dependent function is justified only 



when polarization is neglected (e.g., IFrischl 119961 ). Since resonance polarization is largely 



controlled by the anisotropy of the radiation field, the use of angle-averaged functions may 
not be as suitable for pol arized tran s fer as for non-polarized transfer. Tests of this ap- 
proximation performed by iFaurobertl (119881 ) for purely coherent scattering in semi-infinite 
isothermal atmospheres show that the errors in the linear polarization are very small in the 
wings, while in the line core the polarization peak obtained with the angle- dependent type 
II redistribution function is slightly sharper than with -R^aa- However, such tests need to 
be revised considering more realistic solar model atmospheres. On the other hand, polarized 
radiative transfer computations of the Hanle effect with angle-averag ed PRD have been per- 



formed for isothermal model atmospheres by several researchers (e.g.. lFaurobert-Scholllll991 
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Nagendra et al.l 119991 ; iFluri et al.l 120031 ; iSampoorna et al.l l2008bl ) . Even though the Hanle 
effect is confined to the line core, there are considerable differences in th e linear polariza- 
tion profiles computed usin g CRD and angl e-averaged PRD functi ons (e.g.. iFaurobert-Scholl 
199ll ; iNagendra et al.lll999l ). Furthermore, iNagendra et al.l (120021 . see also Sampoorna et al. 
2008a) have investigated the reliability of angle- averaged functions for the case of Hanle effect 
in isothermal atmospheres, showing that the Stokes U parameter is relatively more sensitive 
than Stokes Q regarding the use of angle-averaged versus angle- dependent functions. There- 
fore, to further assess the relative importance of PRD effects, we need to carry out detailed 
investigations similar to the one reported in this paper but including the Hanle effect of weak 
magnetic fields. Obviously, such radiative transfer problems are significantly more compli- 
ca ted that the one considered he r e, but we think that the same numerical method presented 
in ISampoorna fe Trujillo Buenol ( 120101 ) can be suitably generalized to solve efficiently such 



more general PRD problems, both in ID and 3D geometries. 
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Table 1. The height in the VALC model atmosphere at which lines of different strength r 
have r 0fl = 1 . Other model parameters are A = 5000 A, e = 10" 4 , a = 10~ 3 and 

r E /r R = d^/t r = o. 
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z(r Xfi = 1) in 
x = and /i = 


km 
0.11 


z{t xix = 1) in km 
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Fig. 1. — Temperature stratification of the VALC (dotted line) and MCO (solid line) solar 
model atmospheres. Note that between 600 and 1300 km (hereafter, the "lower chromo- 
sphere") the MCO model is cooler than VALC. 
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Fig. 2. — The continuum optical depth, r c versus height in the VALC model atmosphere 
(solid line) and for a exponentially stratified atmosphere (dotted line). In this paper we have 
used t c = 2.2exp(—z/H), with H = 120km and z in km. 
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Fig. 3. — The anisotropy versus the reduced frequency x at the height in the VALC model 
atmosphere where r XfJi = 1 for a LOS with ji = 0.11. Panels (a) and (b) correspond to PRD 
and CRD respectively. The other model parameters are e = 10" 4 , r = 10~ 5 , a = 10~ 3 and 
Fe/Fr = .D^/Tr = 0. In both panels the solid line is A and the dotted line is A. Notice 
that A CRD and A PRD are similar, while A C rd and A PRD differ greatly, particularly in the 
wings. 
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Fig. 4. — Sensitivity of the PRD solutions in the VALC model to the line strength parameter 
r. The different line types are the following. Solid: r = 10" 6 ; dotted: r = 10" 5 ; dashed: 
r = 10" 4 ; dot-dashed: r = 1CT 3 ; dash-triple-dotted: r = 1CT 2 ; long-dashed lines: r = 0.7. 
Other model parameters are A = 5000 A, e = 10~ 4 , a = 10~ 3 and r E /r R = /T R = 0. 
The top solid line in the Si panel is the Planck function B for the temperature stratification 
of the VALC model. The vertical lines in the same panel show the height at which r 0At = 1 
for a LOS with \x = 0.11. The inset in the —Q/I panel shows the line core region in more 
detail. The symbol 7 cont denotes the intensity at very large distance from line center. 
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Fig. 5. — Same as Figure HJ but for CRD. The line types and the model parameters are 
exactly the same as in Figure HI 
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Fig. 6. — The effect of e on PRD solutions in the VALC model. The different line types are 
the following. Solid: e = 10~ 5 ; dotted: e = 1(T 4 ; dashed: e = 10~ 3 ; dot-dashed: e = 10~ 2 . 
Other model parameters are A = 5000 A, r = 10~ 5 , a = 10" 3 and T E /r R = D<® /F R = 0. 
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Fig. 7. — Same as Figure El but for CRD. The line types and the model parameters are 
exactly the same as in Figure El 
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Fig. 8. — Sensitivity of the PRD solutions in the VALC model to variations in the damping 
parameter a. The different line types are the following. Solid: a = fCT 6 ; dotted: a = 1CT 4 ; 
dashed: a = fCT 3 ; dot-dashed: a = 10~ 2 ; dash-triple-dotted: a = 0.1. Other model 
parameters are A = 5000 A, e = 10" 4 , r = 10" 5 , and r E /r R = £> (2) / r R = 0. 
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Fig. 9. — Same as Figure El but for CRD. The line types and the model parameters are 
exactly the same as in Figure El 
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Fig. 10. — The effect of the elastic collisional rate I^e/I^r on the PRD solutions in the VALC 
model. The different line types are the following. Solid: Te/Tr = 0; dotted: Te/Tr = 0.1; 
dashed: Te/Fr = 0.25; dot-dashed: I^e/Tr = 1; dash-triple-dotted: I^e/Tr = 5; long- 
dashed: Te/Fr = 10. Other model parameters are A = 5000A, e = 10~ 4 , r = 10" 5 , 
a = 10~ 3 and /F R = 0.5r E /r R . 
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Fig. 11. — Same as Figure ITUl but for CRD. The line types and the model parameters are 
exactly the same as in Figure [TUl 
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Fig. 12. — PRD solutions in the VALC model for two spectral lines with different wave- 
lengths. Solid line: A = 3933 A; dotted line: A = 5000 A. Other model parameters are 
e = 10~ 4 , r = 10~ 5 , a = 10" 3 and r E /r R = /T R = 0. 
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Fig. 13. — PRD solutions in the VALC model for spectral lines of the upper chromosphere 
and transition region. The different line types are the following. Solid: r = 10~ 6 ; dotted: 
r = 5 x 10~ 7 ; dashed : r = 4 x 10" 7 ; dot-dashed : r = 3 x 1CT 7 ; dash-triple-dotted : 
r = 2.5 x 1CT 7 ; long-dashed lines: r = 2 x 1CT 7 . Other model parameters are A = 5000 A, 
e = 10" 4 , a = 10" 3 and r E /r R = /T R = 0. 
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Fig. 14. — The effect of the W2 factor on the PRD (top panels) and CRD (bottom panels) 
solutions in the VALC model for a strong line with r = 2 x 10~ 7 . Solid line : W2 = 1 (which 
corresponds to a line transition with J; = and J u = 1); dotted line: W2 = 0.5 (which 
corresponds to a line transition with J\ = 1/2 and J u = 3/2). Other model parameters are 
A = 5000 A, e = 10" 4 , a = 10~ 3 and r E /r R = /T R = 0. 
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Fig. 15. — The effect of the atmospheric thermal structure on the PRD solutions. Solid 
lines : MCO model, dotted lines : VALC model. Other model parameters are A = 5000 A, 
e = 1(T 4 , r = 10" 5 , a = lO" 3 and T E /r R = /T R = 0. 
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Fig. 16. — The effect of the elastic collisional rate 1^/1^ on the PRD solutions in the MCO 
model. The different line types are the following. Solid: Te/Fr = 0; dotted: Te/Fr = 0.1; 
dashed: I^e/Tr = 0.25; dot-dashed: Te/Fr = 1; dash-triple-dotted: Te/Fr = 5; long- 
dashed lines: Te/Fr = 10. Other model parameters are A = 5000 A, e = 10~ 4 , r = 10" 5 , 
a = 1Q- 3 and D^/T K = 0.5r E /r R . 



